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Abstract — We propose a strange-attractor model of tumor 
growth and metastasis. It is a 4-dimensional spatio-temporal 
cancer model with strong nonlinear couplings. Even the 
same type of tumor is different in every patient both in 
size and appearance, as well as in temporal behavior. This 
is clearly a characteristic of dynamical systems sensitive to 
initial conditions. The new chaotic model of tumor growth and 
decay is biologically motivated. It has been developed as a 
live Mathematica demonstration, see Wolfram Demonstrator site: 



If tumor growth is chaotic, this could explain the unreli- 
ability of treatment and prediction of tumor evolution. More 
importantly, if chaos is established, this could be used to adjust 
strategies for fighting cancer using chaos control and/or anti- 
control^ 

In this paper, a plausible chaotic diffusion model of tumor 
growth and decay is presented, as a combination of theoretical 
modelling and empirical search. 



http://demonstrations.wolfram.com/ChaoticAttractorInTumorGrowth/ 
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Index Terms — Reaction-diffusion tumor growth model, chaotic 
attractor, sensitive dependence on initial tumor characteristics 

I. Introduction 

Cancer is one of the main causes of morbidity and mortality 
in the world. There are several different stages in the growth 
of a tumor before it becomes so large that it causes the 
patient to die or reduces permanently their quality of life. 
Developed countries are investing large sums of money into 
cancer research in order to find cures and improve existing 
treatments. In comparison to molecular biology, cell biology, 
and drug delivery research, mathematics has so far contributed 
relatively little to the area QJ. 

A number of mathematical models of avascular (solid) 
tumor growth were reviewed in [2|. These were generally 
divided into continuum cell population models described by 
diffusion partial differential equations (PDEs) of continuum 
mechanics 0, combined with chemical kinetics, and dis- 
crete cell population models described by ordinary differential 
equations (ODEs). 

On the other hand, in many biological systems it is possible 
to empirically demonstrate the presence of attractors that 
operate starting from different initial conditions 0. Some 
of these attractors are points, some are closed curves, while 
the others have non-integer, fractal dimension and are termed 
"strange attractors" [6|. It has been proposed that a prerequisite 
for proper simulating tumor growth by computer is to establish 
whether typical tumor growth patterns are fractal. The fractal 
dimension of tumor outlines was empirically determined using 
the box-counting method Q. In particular, fractal analysis 
of a breast carcinoma was performed using a morphometric 
method, which is the box-counting method applied to the 
mammogram as well as to the histologic section of a breast 
carcinoma [8]. 
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II. A REACTION-DIFFUSION CANCER GROWTH MODEL 

Recently, a multiscale diffusion cancer-invasion model 
(MDCM) was presented in (9), 03, 02, 01, 03, 03, 
lfl5l . 02), 03, which considers cellular and microenviron- 
mental factors simultaneously and interactively. The model 
was classified as hybrid, since a continuum deterministic 
model (based on a system of reaction-diffusion chemotaxis 
equations) controls the chemical and extracellular matrix 
(ECM) kinetics and a discrete cellular automata-like model 
(based on a biased random-walk model) controls the cell 
migration and interaction. The interactions of the tumor 
cells, matrix-metalloproteinases (MMs), matrix-degradative 
enzymes (MDEs) and oxygen are described by the four 
coupled rate PDEs: 
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where n denotes the tumor cell density, / is the MM- 
concentration, m corresponds to the MDE-concentration, 
and c denotes the oxygen concentration. The four variables, 
n, m, /, c, are all functions of the 3-dimensional spatial vari- 
able x and time t. All equations represent diffusion ex- 
cept (H, which shows only temporal evolution of the MM- 
concentration coupled to the MDE-concentration. D n is the 
tumor cell coefficient, D m is the MDE coefficient and D c > 

'The chaotic behavior of a system may be artificially weakened or sup- 
pressed if it is undesirable. This concept is known as control of chaos. The 
first and the most important method of chaos control is the so-called OGY- 
method, developed by 1 19']. However, in recent years, a non-traditional concept 
of anti-control of chaos has emerged. Here, the non-chaotic dynamical system 
is transformed into a chaotic one by small controlled perturbation so that 
useful properties of a chaotic system can be utilized |20|, |21|. 
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is the oxygen diffusion coefficient, while x, M> \ a > 7> 
are positive constants. The other terms respectively denote: 

XV • (riVf)— haptotaxis; 

fiN— production of MDE by tumor cell; 

Am— decay of MDE; 

5mf— degradation of MM by MDE; 

ac— natural decay of oxygen; 

771— oxygen uptake; and 

(3f— production of oxygen by MM. 

Because of its hybrid nature (cells treated as discrete enti- 
ties and microenvironmental parameters treated as continuous 
concentrations), the 4-dimensional (4D) model (Q3-(|4]i can 
be directly linked to experimental measurements of those 
cellular and microenvironmental parameters recognized by 
cancer biologists as important in cancer invasion. Furthermore, 
the fundamental unit of the model is the cell, and the complex 
collective behavior of the tumor emerges as a consequence 
of interactions between factors influencing the life cycle and 
movement of individual cells lEfl. ifTTI. fl2l. |[T5l. fltl. ifTTl. 

In order to use realistic parameter values, the system of rate 
equations <fXJ — <f4j> was non-dimensionalised. The resulting 4D 
scaled system of rate PDEs IfTTI is given by 
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In ifTTIl the values of the non-dimensional parameters were 
given as: 

d n = 0.0005, d m = 0.0005, d c = 0.5, p = 0.01, (9) 
Tj = 50, k = 1, a = 0, v = 0.5, uj = 0.57, c/> = 0.025. 

The 4D hybrid PDE-model ([T])-© can be seen as a special 
case of a general multi-phase tumor growth PDE (j2| equation 
(12)), 



V-(A$i)+A*(*i, Ci)-Hi(*i, Q), (10) 



where for phase i, $i is the volume fraction $i = 1), 
Vj is the velocity, Di is the random motility or diffusion, 
Ai($i,Ci) is the chemical and phase dependent production, 
and p,i($>i, Ci) is the chemical and phase dependent degrada- 
tion/death. The multi-phase model (fTOt has been derived from 
the basic conservation equations for the different chemical 
species, 

i)( 

V-N,- = R. 



dt 



where Ci are the concentrations of the chemical species, 
subindex a for oxygen, 6 for glucose, c for lactate ion, d for 
carbon dioxide, e for bicarbonate ion, / for chloride ion, and 
g for hydrogen ion concentration; Nj is the flux of each of the 
chemical species inside the tumor spheroid; and Pi is the net 
rate of consumption/production of the chemical species both 
by tumor cells and due to the chemical reactions with other 
species. 



In this paper we will search for a temporal 3D cancer 
chaotic attractor 'buried' within the 4D hybrid spatio-temporal 
model ©-©. 

III. A CHAOTIC MULTI-SCALE CANCER GROWTH/DECAY 
MODEL 

From the non-dimensional spatio-temporal AC model (|5]l- 
dHJ, discretization was performed by neglecting all the spatial 
derivatives resulting in the following simple 4D temporal 
dynamical system. 
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When simulated, the temporal system (fTTl i- (fT4l i with the set 
of parameters (O exhibits a virtually linear temporal behavior 
with almost no coupling between the four concentrations that 
have very different quantitative values (all phase plots between 
the four concentrations, not shown here, are virtually one- 
dimensional). To see if a modified version of the system (fTTl l- 
(fl4] i could lead to a chaotic description of tumor growth, four 
new parameters, a, (3, 7, and 5 were introduced. The resulting 
model is 
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The introduction of the parameters (a, j3, 7, S) was moti- 
vated by the fact that tumor cell shape represents a visual 
manifestation of an underlying balance of forces and chem- 
ical reactions [32 1 . Specifically, the parameters represent the 
following quantities: 

a = tumor cell volume (proliferation/non-proliferation 
fraction), 

= glucose level, 

7 = number of tumor cells, 

6 = diffusion from the surface (saturation level). 

A tumor is composed of proliferating (P) and quiescent (or 
non-proliferating) (Q) cells. Tumor cells shift from class P to 
class Q as the tumor grows in size [33). Model dependence 
on the ratio of proliferation to non-proliferation is introduced 
via the first parameter, a. The discretization of equation (0 
leads to cell density being modelled as a constant in equation 
(fTTl i. Accordingly, cell density does not play a role in the 
dynamics. In (fT~5T> — (TTHl the cell density is re-introduced into 
the dynamics via the cell number, 7. The importance of 
introducing 7 also appears in connection with the cyclin- 
dependent kinase (Cdk) inhibitor p27, the level and activity 
of which increase in response to cell density. Levels and 
activity of Cdk inhibitor p27 also increase with differentiation 
following loss of adhesion to the ECM |29|. 

The ability to estimate the growth pattern of an individual 
tumor cell type on the basis of morphological measurements 
should have general applicability in cellular investigations, 
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cell-growth kinetics, cell transformation and morphogenesis 

m. 

Cell spreading alone is conducive to proliferation and in- 
creases in DNA synthesis, indicating that cell morphology is a 
critical determinant of cell function, at least in the presence of 
optimal growth factors and extracellular matrix (ECM) binding 
||35l . In many cells, the changes in morphology can stimulate 
cell proliferation through integrin-mediated signaling, indi- 
cating that cell shape may govern how individual cells will 
respond to chemical signals |36|. 
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Fig. 1, A 3D Lorenz-like chaotic attractor (see [221, 1231) from the modified 
tumor growth model I15M18I . The attractor effectively couples the MM- 
concentration /, the MDE-concentration in, and the oxygen concentration c 
in a mask-like fashion. 

Parameters (a, (3,^,6), introduced in connection with can- 
cer cells morphology and dynamics could also influence the 
very important factor chromatin associated with aggressive 
tumor phenotype and shorter patient survival time. For com- 
putations, the parameters were set to a = 0.06, (3 = 0.05, 7 
= 26.5 and 5 = 40. Small variation of these chosen values 
would not affect the qualitative behavior of the new temporal 
model (fl5]l-(fT8]l. Simulations of <fT3T> — (fl~8T>, using the same 
initial conditions and the same non-dimensional parameters 
as before, show chaotic behavior in the form of Lorenz-like 
strange attractor in the 3D (/ — m — c) subspace of the full 
4D (n — / — m — c) phase-space (see Figure [TJ. 



IV. A CHAOTIC REACTION-DIFFUSION CANCER 
GROWTH/DECAY MODEL 

The spatio-temporal system of rate PDEs corresponding to 
the system in (15)— (18) provides the following multi-scale 



cancer invasion model. 
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The new tumor-growth model (|T9b— (f22b retains all the 
qualities of the original AC model (H)-© plus includes the 
temporal chaotic 'butterfly'-attractor. This chaotic behavior 
may be a more realistic view on the tumor growth, including 
stochastic-like long-term unpredictability and uncontrollabil- 
ity, as well as sensitive dependence of a tumor growth on its 
initial conditions. 

Based on the new tumor growth model (fT9l-(l22l. we have 
re-formulated the following generic concept of carcinogenesis 
and metastasis (see scheme on the following page). With the 
model of strongly coupled PDEs, remodelling of extracellular 
matrix (ECM) causes a whole process of the movement of 
invading cells with increased haptotaxis, and at the same time 
decreased enzyme productions level. This can alter chromatin 
structure, which plays an important role in initiating, propa- 
gating and terminating cellular response to DNA damage [28]. 
The effect of haptotaxis in the process of cells invading could 
be modelled with travelling-wave (Fisher) equation [18|. 

The proposed model (fl~9b — (f22l> describes chaotic behavior 
relevant to the invasion of cancer cells. As devices for con- 
trolling metastasis/chaos we suggest the following processes: 
Cellular retraining of cancer stem cells and/or activation of 
positive function of cyclin-dependent kinase inhibitor p27 
and/or decreased expression of SATB1, which is correlated 
with aggressive tumor phenotype in breast cancer and shorter 
patient survival time ll27l . (29 1, [30|, [31 1. To connect the new 
parameters (a, f3, 7, S) to the therapeutic regime we note the 
fact that significant factor of any therapy is tumor re-growth 
during the rest periods between therapy applications, which is 
again dependent on proliferation fraction dynamics a ll33l . 

Also, the age of the tumor may translate to different levels 
of response to a standard therapy. As mentioned earlier, the 
number of cells, 7, is connected with the Cdk inhibitor p27, 
for which levels and activity increase in response with the 
cell density 7, differentiation following loss of adhesion to 
the ECM. Cdk inhibitor p27 regulates cell proliferation, cell 
motility and apoptosis and is the essential element for under- 
standing transduction pathways in the regulation of normal 
and malignant cell proliferation as well as it is new hope for 
therapeutic intervention [29 1 . 

Introducing the parameter set (a,/3,7, 6) into the model 
of tumor cell morphology and function could lead to insight 
into the relationship between these parameters and chromatin. 
By varying these parameters, it may be possible to predict 
situations that result in dangerous levels of chromatin. High 
levels of chromatin are associated with cancer metastasis and 
some of the most aggressive cancer types l30ll . (3T|. 

The new temporal model (fT~5T> — (fT~8~b (as well as the corre- 
sponding spatio-temporal tumor-growth model (fT9l-(f22l) is 
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not sensitive to variation of the (a, f3, 7, 6)— values, but is 
sensitive to their corresponding places in the equations. 

V. Conclusion 

A plausible chaotic multi-scale cancer-invasion and decay 
model has been presented. The new model was formulated by 
introducing nonlinear coupling into the existing hybrid mul- 
tiscale Anderson-Chaplain model. The new model describes 
chaotic behavior, as well as sensitive dependence of a tumor 
evolution on its initial conditions. Effective cancer control is 
reflected in progressive reduction in cancer mortality [27|. The 
proposed model suggests a possible solution to carcinogenesis 
and metastasis, by combining mathematical modelling with 
latest medical discoveries. 
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